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I.  IN.tfiQPyCII.ON 

Velocity  sheared  nixing  layers,  or  shear  layers,  are 
en  important  design  consideration  for  many  high  power 
lasers.  Mixing  layers  will  occur  across  openings  in  the 
skin  of  elrcraft,  for  example,  where  some  lasers  may  have 
application.  And,  they  are  found  inside  the  laser  itself 
in  many  instances:  they  can  be  generated  in  the  laser 
resonator,  and  they  are  significant  in  aerodynamic 
windows^ . 

In  addition,  the  mixing  layer  is  Interesting  in  itself 
as  the  subject  for  reaeerch  in  fluid  dynamics.  Many  recent 
research  papers  have  examined  the  effects  on  light  as  it 
passes  through  a  mixing  layer.  In  some  cases,  a  laser  beam 
is  used  as  a  diagnostic  tool  to  determine  mixing  layer 
characteristics  based  on  the  distortion  in  the  near  field 
of  the  beam  after  passing  through  the  layer.  In  other 
cases  it  is  the  distortion  itself  which  is  the  subject  of 
research.  Poling^  (1985),  for  example,  used  near  field 
Interferometry  (among  other  techniques)  to  study  mixing 
layer  spreading  rates  and  density  gredients. 
Christiansen  et  al^  (1984)  did  similar  work,  but  also 
determined  the  mean  square  phase  error  for  use  in 
calculating  the  Strehl  ratio4,  I/I0»  which  is  a  useful 
quantity  for  measuring  the  degradation  of  the  laser  beam. 


Legner  at  al5  (1978)  and  Vu  at  al8  (1980)  developed  models 
for  tha  raductlon  in  beam  intensity  at  the  peak  of  the  far 
field  distribution  and  obtained  experimental  data  for  the 
velocity  matched  mixing  layer  similar  to  that  found  in  the 
aerodynamic  window  of  a  CO2  laser.  Wasserstrom  at  al? 
demonstrated  theoretically  that  a  ay stem  of  crossed  jets 
(and  their  attendant  shear  layers)  could  be  used  to 
partially  compensate  for  phase  error  caused  by  thermal 
blooming.  In  each  of  these,  interf erograms  of  the  near 
field  are  the  primary  source  of  data  on  the  distortion  of 
the  beam.  From  near  field  data  one  may  compute  the 
degradation  that  would  be  seen  in  the  far  field. 

Geargeoura8  (1980),  on  the  other  hand,  measured  the 
the  far  field  degradation  directly.  By  inserting  various 
size  circular  apertures  into  the  beam  after  it  had  passed 
through  the  mixing  layer,  he  was  able  to  measure  the  the 
average  cross  sectional  distribution  of  beam  power  and 
compare  it  to  the  undistorted  case. 

This  thesis  describes  a  different  technique  that  was 
used  to  obtain  data  directly  from  the  far  field  of  a  laaer 
beam  after  paaaing  through  a  mixing  layer.  Further,  it 
presents  dats  obtained  by  the  technique  on  a  variety  of 
mixing  layers  at  two  downatraam  locations  and  using  two 
beam  diameters.  The  work  presented  here  is  very  similar  to 
that  of  Geargeoura  in  that  tha  far  field  intensity 
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distribution  is  measured  directly  for  circular  beams 
passing  through  mixing  layers  generated  by  jets  of  various 
gases.  The  primary  difference  being  the  technique  for 
measuring  the  far  field  intensity  distribution  of  the  bean. 
Geargeoura  used  a  mask  to  sample  various  portions  of  the 
degraded  beam  and  a  lens  to  focus  the  sampled  portion  onto 
a  single  photodiode  to  measure  the  power  in  that  portion. 
Here,  the  majority  of  the  beam's  far  field  is  captured 
simultaneously  by  a  charge  coupled  device  (CCD>  camera 
array  of  256  x  256  photodiodes.  Other  than  the 
demagnif icatlon  of  the  beam  to  shorten  the  distance  to  the 
far  field,  similar  to  Geargeoura's,  no  lens  is  used  to 
focus  the  light  being  sampled.  This  technique  does  not 
suffer  from  the  requirement  for  a  separate  run  for  each 
region  of  the  beam  to  be  examined  since  it  samples  almost 
the  entire  beam  at  once.  Nor  does  it  depend  on  the  mask 
being  centered  with  respect  to  the  distorted  beam  since  it 
does  not  use  this  type  of  mask. 

This  thesis  covers,  in  detail,  the  procedures  used  to 
obtain  far  field  diffraction  pattern  images  and  their 
subsequent  analysis.  Even  though  the  conditions  addressed 
here  were  limited  to  zero  velocity  ratio  layers  with  a 
range  of  density  ratios  from  0.195  to  7.22  at  a  Mach  number 
of  0.1,  it  is  the  first  step  in  developing  this  technique 
to  cover  the  wider  range  of  conditions  of  interest. 
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surrounding  gas,  ths  boundary  batwaen  the  two  gasas  ia 


unstabla  bscauss  of  ths  vslocity  discontinuity  and 


rssulting  shsar  atrssa  at  ths  boundary.  Vary  quickly,  ths 


boundary  "rolls  up"  into  vorticsa  which  six  ths  various 


propsrtiss  of  the  two  gasas.  As  the  vortices  progress  down 


stream  they  pair  and  combine  with  other  vortices  end  grow 


in  size,  thus  the  region  of  mixing  grows  in  thickness10. 


This  growing  region  of  vortices  moving  down  stream  la  the 


turbulent,  velocity  sheared  mixing  layer.  A  sketch  of  auch 


a  layer  is  shown  in  figure  1. 


The  distance  between  the  vortices  as  they  ere  formed. 


and  thus  an  indication  of  how  far  downstream  from  the 


nozzle  they  form,  is  expected  to  be  proportioal  to  boundary 


layer  thickness  at  the  nozzle  exit11,  and  by  lOOO  momentum 


thicknesses  downstream  the  flow  is  self  similar12. 


Previous  researchers2  using  the  same  nozzle  that  is  used 


here  have  calculated  the  momentum  thlcknesa  of  the  boundary 


layer.  Using  the  Pohlhausen  analysis  with  Walz'a  solution. 


as  discussed  in  Schlichting1^,  for  Hach  number  M  *  0.1  the 


momentum  thickness  of  the  boundary  layer  at  the  nozzle  exit 


has  been  celculated  to  be  0.0028  cm.  Thus,  the  region  of 


vortices  should  be  well  established  1  or  2  cm  downstream. 


VJV.V 
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Fig.  i.  Mixing  Layer.  Sketch  of  a  aixing  layer 
caueed  by  a  free  jet. 


7 


i 

i 

One  of  the  properties  mixed  by  the  vortices  is  the 
index  of  refraction.  For  a  gas,  the  index  of  refraction 
may  be  obtained  from14: 

P 

n  =  1  +  P  — 

ps  (1) 

Where 

P  =  a  tabulated  constant 
p  =  density 

Ps  *  density  at  standard  conditions 
n  *  index  of  refraction. 


and  the  tabulated  constant,  15 ,  depends  on  the  gas  and  the 
wavelength (a)  of  light  being  used. 

Brown  and  Roshko10  demonstrated  that,  as  predicted  by 
theory13,  the  thickness  of  the  mixing  layer  grows  linearly 
with  distance  downstream,  hence  the  mixing  layer  shape  may 
be  approximated  as  a  wedge.  According  to  Vu  et.  al.^  the 
thickness  of  this  wedge,  L,  may  be  expressed  as: 


L  =  c 


’u.-iV 

U.  +  U,, 


<  2) 


where: 


x  *  distance  downstream. 

Ul  *  x-velocity  of  high  speed  stream. 

U2  *  x-velocity  of  low  speed  stream  (zero  in  this 

case) . 

c  »  a  constant  for  each  nixing  layer. 

Brown  and  Foshko^O  measured  c  and  found  it  to  be  a 
weak  function  of  density  ratio  having  values  of  0.51,  0.38, 
and  0.28  corresponding  to  density  ratios  of  P»/ P i  *  7,  1, 

and  1/7  respectively. 

It  is  known  however,  that  for  the  nozzle  used  here, 
the  linear  growth  approximation  given  by  eqn.  (2)  is  only 
valid  up  to  a  distance  of  about  2  exit  widths,  or  about 
2.8  cm  downstream2.  Nevertheless,  if  the  layer  is 

represented  as  a  wedge  shaped  region  having  an  index  of 
refraction  equal  to  the  mean  of  the  index  of  refraction  of 
the  two  gases  (air  and  the  jet>,  then  the  beam  would  be 
deflected  according  to  Snell's  law  at  each  of  the  two 
interfaces  (see  fig.  2): 


sin  (  sin  (  0t ) 


(3) 


Fig.  2.  Deflection  of  e  light  ray  passing  through  a 
wedge  aheped  region  whose  refractive  index  is  the 
average  of  the  refractive  indices  of  the  regions  on 
either  side  of  the  wedge. 
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Thus,  the  total  deflection  of  the  beam  is  given  by 


&  =  J  arctan 


'  dL  ' 


l  *c  J 


-  0 


t  2 


Where: 


0t2=  aroin 


\  n£H 

1  2  rij 


sln(y-0t  j) 


) 


j  =arcsin 


r  2na 

1  VnJ 


y 

—  +  ♦ 
2 


el2  =  y-0 


t  i 


and. 


Y  *  arctan  (dL/dx),  the  wedge  angle. 
na  3  refractive  index  of  air. 
nj  3  refractive  index  of  the  jet. 

6  =  angular  deflection  of  the  ray. 

<p  3  “tilt"  of  the  wedge  toward  the  incoming  beam 
(aae  fig.  2> 

The  values  of  6  predicted  by  eqn.  (4)  will  be  mentioned 
briefly  in  Diacuaaion  and  Conclusion  sections. 
Surprisingly  enough  though,  eqn.  (4>  predicts  that,  for  a 
homogeneous  wedge,  when  $  is  zero  the  angular  deflection  of 
the  beam  is  essentially  zero  as  well. 

At  best  however,  eqn.  (4>  is  an  estimate  of  the  time 
averaged  deflection  of  the  beam.  Because  the  vortices  in 
the  mixing  layer  are  traveling  through  the  beam  as  they 
move  downstream  they  cause  nonsteady  effects  as  well.  The 
vortices  move  at  a  speed  which  is  the  average  between  that 
of  the  jet  and  that  of  the  air.  Since  the  air  is 

stationary,  the  vortices  generated  by  a  jet  with  N  ■  0.1, 
for  example,  move  downstream  at  roughly  1700  cm/sec  if  the 
jet  gas  is  nitrogen.  For  a  helium  jet  their  speed  would  be 
about  5000  cm/sec,  and  roughly  1300  cm/sec  for  a  carbon 
dioxide  jet. 

In  addition,  the  actual  density  distribution  in  the 
vortices  is  highly  nonlinear.  This  means  that  the  density 


distribution  across  tha  mixing  layer  is  nonlinear.  This 
causes  wavefront  distortion  to  be  added  to  the  beam 
deflection  or  tilt.  Vu  at.  al.&  have  derived  the  following 
theoretical  relationship  to  predict  this  degradation  of  the 

beam: 


Io-I 

I 


=  J  kJ  An*  ar1  c* 


Ut  -  u» 

.  U,  +  Uj  J 


(5) 


In  this  expression: 

I0  «  nondegraded  beam  intensity  at  the 
center  of  the  far  field  focus. 

I  *  degraded  beam  intensity  at  the  center  of 
the  far  field  focus. 

k  *  wave  number 

An  *  index  of  refraction  difference  across 
the  mixing  layer. 

o  ■  a  measured  quantity. 

The  quantity  a  is  proportional  to  the  turbulence 


intensity  and  has  a  value  close  to  that  of  the  density  or 
temperature  fluctuations: 
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Batt1^  has  measured  ths  temperature  fluctuation  to  bs 
approximately  0.15,  while  Fiedler^^  measured  it  to  be 
nearly  0.20. 

In  his  1969  paper,  Sutton*?  presented  a  relationship 
related  to  that  in  eqn.  (5).  In  fact,  the  term  on  the 
right  in  eqn.  (5)  may  be  regarded  as  the  first  term  in  s 
Taylor  series  expansion  of  the  right  side  of  Sutton's  1969 
relationship: 


Io-I 

=  1  -  exp 

-  *  k*  An*  a*  c* 

C 

i 

c 

N 

_ > 

2  ■ 

X* 

Io 

„ 

l  U,  +  U2  J 

« 

The  left  side  of  eqns.  (5)  and  (6>  is  just  one  minus 
the  Strehl  ratio.  The  Strehl  ratios  obtained  from  eqn.  (6) 
will  be  compared  to  measurements  and  discussed  briefly  in 
the  Discussion  and  Conclusions  sections. 


Ill 


In  order  to  measure  the  beam  degradation  for  thia 
experiment,  the  equipment  arrangement  sketched  in  figure  3 
was  used.  The  beam  waa  spatially  filtered,  expanded, 
collimated,  and  then  masked  by  a  circular  aperture  to 
produce  a  parallel  beam  of  nearly  constant  intensity 
distribution  across  a  circular  cross  section.  This  beam 
was  then  passed  through  the  mixing  layer  of  a  jet  exiting  a 
rectangular  nozzle  into  a  three  sided  rectangular  test 
section  with  optical  quality  glass  walla.  The  beam  then 
entered  a  reducing  telescope  of  sufficient  power  to  bring 
the  Fraunhoffer  far  field  diffraction  pattern  within 
convenient  range  of  the  optical  bench.  A  charge  coupled 
device  (CCD)  electronic  camera  then  captured  the  image  of 
the  diffraction  pattern  for  digital  storage  and  analysis  on 
a  dedicated  microcomputer.  A  detailed  description  of  the 
equipment  and  experimental  setup  is  given  below. 


>nd  Teat  Section; 


A  custom  built  nozzle  was  used  to  produce  the  mixing 
layer.  The  nozzle  waa  operated  at  a  nominal  Mach  number  of 
0.1  with  a  variety  of  gases  to  generate  mixing  layers,  in 
room  air,  having  density  ratios  from  .195  to  7.22.  The 
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Fig.  3.  Sketch  of  the  equipment  arrangement.  ET  is 
en  expending  telescope  with  epetiel  filter.  M  la  the 
aperture  aaak.  TS  la  the  test  section.  RT  is  the 
reducing  telescope  without  spatial  filter.  A  is  the 
stack  of  attenuating  neutral  density  filters.  BG  la 
the  infrared  filter.  C  is  the  CCD  camera.  V  la  the 
MicroVAX  computer  with  Microtex  camera  driver  and  A/D 
board . 


uuuA«nuiRfiuMuniimm 
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nozzle  exit  cross  section  was  rectangular,  1.4  cm  by 
1.4  cm,  and  exited  into  a  rectangular  test  section 
consisting  of  optical  quality  glass  walls  on  three  sides 
with  the  fourth  side  open  to  room  sir  (see  fig.  4).  This 
configuration  produced  a  free  gaa/air  mixing  layer  at  the 
open  side  of  the  test  section.  The  laser  beam  entered  the 
test  section  through  the  open  side  and  mixing  layer,  and 
exited  through  the  boundary  layer  and  glass  wall  at  the 
opposite  side.  The  operation  of  the  nozzle  and  quality  of 
the  mixing  layer  were  verified  by  Schlieren  photography. 
For  a  more  detailed  discussion  of  the  nozzle,  see  Poling^. 

Optical  Bench: 

Laser:  The  laser  used  for  this  research  was  an 
unpolarized,  continuous  laser  with  a  rated  output  of  0.9S 
milliwatt . 

Bean  Expander:  The  laser  output  been  was  expsnded 
from  1  mm  to  5  cm  in  diameter  and  collimated  by  beam 
expander  with  a  25  micron  spatial  filter. 


Mask:  before  entering  the  test  section,  the  parallel 
beam  from  the  expander  was  masked  by  a  circular  aperture 
machined  in  a  0.126  cm  piece  of  aluminum. 
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Fig.  4.  Nozzle  and  teat  aectlon.  I  la  the  gaa  Inlet. 
P  la  the  plenum.  SC  la  the  aide  cover.  OP  are 
orifice  platea.  PT  la  a  preaaure  tap.  TS  la  the  teat 
aectlon. 
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Reducing  Telescope:  The  beam  which  exited  the  teat 
section  wee  damagnlfled  by  e  40  power  bees  expender  without 
spatiel  filter  oriented  as  a  reducing  telescope.  The 
telescope  was  focused  to  produce  the  sharpest  possible 
Airy1®  pattern. 

Camera:  The  far  field  pattern  was  placed  directly  on 
the  CCD  array  of  an  EG&G  Reticon  camera  model  9256-B.  This 
camera  has  s  256x256  photodiode  array.  Depending  upon  the 
size  of  the  mask  in  use  and  beam  intensity  at  the 
photodiode  array,  the  array  was  located  between  30  cm  and 
75  cm  from  the  output  lens  of  the  reducing  telescope. 

Optical  Rail:  A  single  extrusion  of  aluminum  channel, 
cut  into  two  pieces,  and  locally  manufactured  carriers, 
which  could  be  locked  to  the  channel  at  any  location  along 
its  length,  were  used  to  mount  the  optical  equipment  on  the 
bench.  One  piece  of  channel,  with  mounts,  supported  the 
laser,  beam  expander,  and  mask  on  the  open  (mixing  layer) 
aide  of  the  test  section.  The  reducing  telescope  and 
camera  were  aimilerly  mounted  on  the  opposite  side  of  the 
test  section.  This  arrangement  placed  all  optical  elements 
in  a  straight  line,  eliminating  the  need  for  beam  steering 
mirrors  and  minimizing  the  distance  required  between 
elements.  Beam  positioning  and  aiming  were  accomplished  by 
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measuring  the  baa*' a  entrance  and  exit  locationa  on  the 
teat  section  end  adjusting  the  input  rail  as  required.  The 
optical  axis  of  the  output  rail  was  then  adjusted  to 
coincide  with  the  axis  of  the  input  rail  already 
eatabliahed. 


Filters:  A  BG-16  filter  was  placed  at  the  camera  to 
reduce  infrared  effects  on  the  recorded  l*age.  In  addition 
various  neutral  density  filters  were  placed  at  the  output 
lens  of  the  reducing  telescope  to  control  the  intensity  of 
the  laser  bee*  before  it  arrived  at  the  photodiode  array. 


Bias  Illumination:  It  was  found  to  be  neceasary  to 
slightly  bias  the  photodiode  array  with  an  additional 
"white"  light  source.  This  was  done  by  placing  a  12  volt 
DC  incandescent  light  source  In  the  field  of  view  of  the 
photodiode  array.  The  illumination  from  this  source  passed 
through  none  of  the  optics  other  than  the  BG-18  infrared 
filter.  This  bias  light  aource  was  powered  by  a  voltage 
regulated  DC  power  supply.  The  distance  from  the  array  to 
this  light  varied  but  was  always  greater  then  1.5  meters. 
The  exact  distance  was  determined  for  each  run  by  the 
combination  of  neutral  density  filters,  array  to  reducing 
telescope  distance,  end  bias  light  location  that  produced  e 
well  aimed  Airy  pattern  with  a  detectable  first  bright  ring 
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and  a  peak  intensity  near  but  not  exceeding  the  saturation 
level  of  the  Microtex  digitizer  (see  Data  Acquisition 
Equipment) .  A  "detectable"  first  bright  ring  was  one  which 
was  apparent  in  the  analog  output  of  the  camera  as  seen  on 
an  oscilloscope  used  for  this  purpose. 


Numerical  data  was  obtained  for  later  analysis  through 
the  use  of  a  Microtex  7405  image  acquisition  system 
consisting  of  an  analog  to  digital  (A/D)  conversion  board 
and  a  logic  board  installed  in  a  Digital  Equipment 
Corporation  MlcroVAX  I.  The  logic  board  provided  all 
necessary  signal  and  power  output  to  operate  the  camera  as 
directed  by  Microtex  provided  FORTRAN  programs  running  in 
the  MlcroVAX  computer.  The  digitized  images  output  by  the 
A/D  converter  were  stored  as  data  files  in  the  MicroVMS 
operating  system  on  the  computer's  herd  disk.  These  data 
files  contained  integer  values  between  0  end  4096, 
representing  the  light  intensity  measured  by  each  pixel  in 
the  256x256  array. 


Four  aeparate  checka  war#  uaad  to  verify 
that  tha  optica  uaad  wara  diffraction  liaitad  and  that  tha 
image  collactad  by  tha  CCD  array  waa  lndaad  tha  Airy 
pattarn  of  the  diffraction  aaak. 

Firat,  ainca  tha  largaat  diffraction  aaak 
aperture  uaad  waa  1  cm  in  diameter,  additional  diffraction 
would  not  be  induced  by  the  4.9  cm  diameter  large  lena  of 
the  reducing  teleacope.  And,  with  tha  objective  lana  baing 
uaad  the  teleacope  haa  a  power  of  40,  thua  the  waiat  of  the 
reduced  beam  may  be  eatimated  aa  0.25  mm,  located  at  the 
email  objective  lena  when  the  teleacope  ia  focuaed.  Since 
thla  email  lena  haa  a  diameter  of  2.5  mm,  it  ia  10  tiaea 
tha  diameter  of  the  beam  waiat  and  it  adda  no  diffraction 
to  the  beam  either. 

Second,  the  diatance,  0,  from  the  output  of 
tha  reducing  teleacope  to  the  far  field  waa  eatimated 
uaing19: 
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where: 

X  *  wavelength, 
a  *  radius  of  tha  baam  waiat. 

D  ■  diatanca  froa  tha  raducing  talascope  to 
tha  film. 

Sinca  tha  diatanca  froa  tha  talaacopa  output  waa  always 
greatar  than  30  cm  for  0.S  ca  diameter  apertures  and  50  cm 
for  1  ca  apertures,  tha  caaara  waa  safely  in  the  far  field 
if  the  reducing  telescope  was  properly  focused. 

Third,  a  photographic  iaage  of  tha  focused  output  of 
the  teleacope  waa  analyzed  to  determine  that  the  locations 
of  tha  relative  alniaa  and  maxima  correspond  to  their 
theoretical  locations  for  Airy  patterns.  Figure  5  is  such 
a  photograph.  The  aperture  diameter  was  0.5  cm  and  the 
distance  of  the  film  from  the  output  lens  of  tha  40  power 
reducing  telescope  waa  201  cm.  Tha  measured  radius  to  the 
center  of  the  first  dark  ring,  ri,  is  approximately 
1.25  cm.  A  radius  of  1.24  cm  ia  predicted  by  the 

equation16: 

DX 

r«  =  1.22  (8> 

2a 

Fourth,  numerical  integration  of  tha  intensity 
distribution  in  recorded  images,  as  a  function  of  radius. 


Fig.  5.  Airy  pattern  for  a  0.5  cm  diameter  aperture 
et  a  distance  of  201  cm  from  the  output  of  a  40  power 
reducing  telescope. 


la* 


24 

was  compared  to  the  theoretical  prediction*®  of  the 
Integration  of  the  intenaity  pattern  of  Airy  patterna. 
Figures  6  and  7a  represent  the  “beat"  and  "worst"  of  the 
reaulta  of  such  comparisons,  respectively.  In  figures  6 
and  7  the  radius  is  represented  by  the  nondimensional 
radius  x  *  kaw,  where  k  la  the  wave  number,  a  is  the  radius 
of  the  beam  waist,  and  w  is  the  ratio  of  the  actual  radius 
and  the  distance  from  the  reducing  telescope. 

Because  the  numerical  integration  of  the  camera  pixel 
response  was,  by  necessity,  normalized  to  the  total  of  the 
responses  of  pixels  which  fell  within  the  largest  circle 
that  would  fit  on  the  array,  plots  like  the  one  in  figure 
7a  were  very  common.  Notice  however,  that  figure  6 
encompasses  three  dark  rings  as  Indicated  by  the  three  flat 
spots  in  the  theoretical  curves.  This  indicates  that  the 
image  in  figure  6  has  received  93. 8X  of  the  total  energy  in 
the  original  beam*®.  Figure  7a,  on  the  other  hand, 
encompasses  only  two  dark  rings  indicating  it  has  received 
91k  of  the  original  beam  energy.  Figure  7b  shows  the 
comparison  of  theory  to  the  numerical  integration  of  the 
data  when  the  numerical  intagration  is  normalized  to  the 
sum  of  the  pixel  responses  divided  by  0.91,  the  result  is 
much  nearer  the  theoretical  curve.  This  further  supports 
the  assertion  that  the  optica  were  diffraction  limited. 
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Fig.  6.  Intensity  integration.  Plot  of  the  numerical 
integration  of  the  pixel  intensity  response  to  an 
undistorted  Airy  pattern  and  theory.  The  theoretical 
curve  rises  faster  at  first,  but  reaches  a  lower 
maximum  than  the  numerical  integration  of  the  data. 


Fig.  7.  Normalizing  effect  on  intensity  integration. 
Plot  of  the  numerical  integration  of  the  pixel 
intensity  response  to  an  undiatorted  Airy  pattern  and 
theory.  The  data  in  7a  and  7b  have  been  normalized  by 
different  values. 

7a:  The  numerical  integration  of  actual  data  rises  to 
a  higher  value  than  the  theory  and  approaches  a  value 
of  1.0,  because  it  is  normalized  to  the  total  light 
actually  detected  by  the  diode  array.  The  theoretical 
curve  however,  la  normalized  to  the  total  light 
contained  in  the  beam. 


7b:  When  the  numerical  Integration  of  the  actual  data 
is  multiplied  by  the  theoretical  ratio  of  the  light 
contained  within  the  second  dark  ring  to  the  total 
light  in  the  beam,  the  result  much  more  closely 


resembles  the  theoretical  curve 


Early  experiment*  with  the  camera  however,  showed  that 
while  the  array  reaponae  waa  fairly  linear  aa  a  whole,  the 
level  of  dark  current  and  random  nolae  waa  not  small  in 
relation  to  the  low  signal  strengths  of  the  subsidiary 
maxima  of  the  Airy  patterns.  Worse,  the  output  of  the 
array  was  not  "flat"  across  its  spatial  dimensions.  That 
is,  the  output  of  array  elements  or  pixels  near  the 
horizontal  center  of  the  array  had  lower  dark  output  levels 
than  did  pixels  near  either  side  of  the  array,  and  there 
was  a  definite  high-low  variation  between  even  and  odd 
numbered  lines  of  pixels  from  top  to  bottom  of  the  array. 
Finally,  even  though  moat  of  the  output  reaponae  of  each  of 
the  pixels  was  fairly  linear,  at  low  input  light  levels 
thia  was  not  ao  and  waa  not  the  same  for  each  pixel. 
Because  of  these  problems,  a  diffuse  "white"  light  source 
was  added  to  the  optical  table.  It  consisted  of  a  pair  of 
6  volt  Incandescent  bulbs  powered  in  series  by  a  voltage 
regulated,  12  volt  DC  power  supply.  Thia  additional  light 
had  the  two  beneficial  effects  of:  raising  the  reference 
output  level  of  all  pixels  with  no  laser  beam  input  to  a 
nearly  equal  level  acroaa  the  array,  and  when  added  to  the 
diffraction  pattern  of  the  laser  beam,  raising  the  first 
subsidiary  maximum  of  the  Airy  pattern  above  the  level  of 
the  dark  current  and  random  noise.  The  intensity  of  the 
additional  light  source  seen  by  the  camera  was  controlled 


by  varying  the  distance  from  the  source  to  the  camera  so 
that  the  output  of  the  camera  with  the  laser  turned  off  was 
nearly  flat  and  with  the  laser  turned  on  did  not  saturate 
the  analog  to  digital  converter  board.  The  distance  froa 
the  bulbs  to  the  camera  was  not  measured  for  each  run  but 
was  on  the  order  of  two  meters,  and  the  maximum  linear 
dimension  of  the  two  bulbs  together  was  approximately  5  cm. 
This  bias  illumination  was  later  subtracted  numerically 
from  the  data.  The  plots  in  figures  &  and  7  were  drawn 
after  the  bias  was  subtracted. 

Acoustic  noise  generated  by  the  operation  of  the 
nozzle  was  also  a  concern.  Previous  interferometry  had 
been  rendered  useless  by  noise  Induced  vibration  in  pieces 
of  optical  equipment.  None  of  the  equipment  that  had 
suffered  from  this  in  the  past  was  used  in  this  experiment, 
but  to  verify  that  the  image  distortion  was  not  vibration 
induced,  two  photographs  were  taken  with  the  beam  passing 
just  outside  the  glass  walls  of  the  test  section.  One 
photograph  was  taken  with  the  jet  turned  off,  the  other 
with  a  helium  jet  running  at  a  Mach  number  of  nearly  one, 
the  most  severe  noise  condition  for  this  nozzle.  These  two 
photographs  are  shown  in  figure  8a  and  b. 


Fig.  8.  The  affect  of  jet  generated  noiae.  Airy 
pattern*  of  a  baaa  paaaing  nearby,  but  not  through  the 
teat  aaction 


8a:  Jet  off  (no  noiae) 


8b:  Helium  jet  on.  Mach  number  nearly  1 
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The  optical  banch  waa  aligned  aa  follows: 

A  sat  of  cross  hairs  wss  scribed  on  an  aluminum  target 
mounted  on  a  slide  mount  that  moved  only  parallel  to  the 
rail  axis.  After  the  laser  was  centered  as  closely  as 
possible  on  the  rail,  aiming  was  done  by  rotating  the  laser 
until  the  unexpended  beam  spot  (approx  1  mm  dla.>  location 
on  the  target  cross  hairs  moved  less  than  one  beam  radius 
when  the  target  waa  moved  from  one  end  of  the  rail  to  the 
other . 

The  beam  expander  wae  then  added  to  the  beam  path. 
The  beam  expander  waa  approximately  aligned  using  back 
reflection  from  a  reflective  cap  fitted  over  the  input 
lena.  Then,  the  input  lena  cap  was  removed  and  a  cap  with 
a  1  mm  diameter  hole  in  it's  center  waa  placed  over  the 
output  lena.  The  expander  was  then  rotated  while  the 
target  waa  slid  back  and  forth  again  until  the  spot  moved 
leas  than  one  spot  radius.  During  this  process,  the  beam 
expander  waa  translated  vertically  or  horizontally,  as 
required,  to  maintain  maximum  intensity  output  aa  estimated 
by  eye.  The  beam  expander  assembly  contained  a  25  micron 
spatial  filter  which  was  adjusted  aa  required  to  maximize 
and  center  its  output  on  the  output  lena,  as  estimated  by 
eye. 
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With  this  done,  all  elements  on  the  rail  were  locked 
in  position  and  the  horizontal  location  of  the  entrance  of 
the  1  mm  diameter  beam  into  the  teat  section  was  measured 
using  a  transparent  plastic  rule.  Then  one  end  of  the  rail 
was  moved  while  the  other  was  held  fixed  so  that  this 
location  was  centered  on  one  side  of  the  test  section. 
That  end  of  the  rail  was  then  fixed  and  the  other  end  moved 
until  the  location  of  the  beam  exit  was  centered  on  the 
opposite  side  of  the  teat  section.  By  repeating  this 
technique  iteratively,  the  beam  path  through  the  test 
section  was  centered. 

The  rail  on  the  output  side  of  the  test  section  was 
then  fastened  to  the  table  approximately  parallel  to  the 
first  rail  by  the  use  of  a  straight  edge.  Final  alignment 
of  the  second  rail  was  done  using  the  sliding  target,  and 
moving  alternate  ends  of  the  rail,  as  before.  The  range  of 
movement  for  the  sliding  target  was  approximately  40  cm  on 
either  rail. 

The  laser  had  previously  been  positioned  so  that  the 
beam  was  normal  to  the  jet  and  approximately  1  cm 
downstream  of  the  nozzle  exit.  Now,  the  beam  expander 
output  lens  cap  was  removed  and  a  circular  diffraction  mask 
with  a  1  mm  diameter  was  inserted  in  the  beam  path  and  its 
vertical  position  adjusted  so  that  the  beam  was  measured  to 
be  1  cm  downstream  of  the  nozzle  exit. 


Next,  the  reducing  teleecope  wee  added  to  the  beam 
path  on  the  opposite  aide  of  the  teat  section  from  the 

laser  and  was  aligned  using  a  lens  cap  having  a  1  am 

diameter  hole  in  its  center  covered  by  a  beam  splitter. 
The  reducing  telescope  was  rotated  and  translated  until  the 
transmitted  portion  of  the  1  mm  beam  passed  through  the 

1  mm  hole  in  the  center  of  the  lenscap  and  the  reflected 

portion  returned  to  the  1  mm  hole  in  the  diffraction  mask. 

Finally,  the  diffraction  mask  mount  was  so  fabricated 
that  masks  with  different  diameters  and  center  locations 
could  be  exchanged  without  disturbing  the  alignment  of  the 
rest  of  the  beam.  With  the  proper  size  and  center  location 
mask  installed,  the  reducing  telescope  was  focused  and  the 
CCD  camera  was  positioned  so  that  an  appropriately  sized 
image  was  centered  on  the  array. 

In  order  to  reduce  the  effects  of  dark  current,  the 
camera  was  actively  cooled  by  attaching  a  heat  sink  to  the 
camera  case  and  forcing  cooled,  dry  nitrogen  through  fins 
in  the  heat  sink.  The  nitrogen  gas  was  cooled  by  passing 
*■  through  s  tube  Immersed  in  liquid  nitrogen.  By  this 
means,  the  camera  case  temperature  was  kept  between  13  and 
15  degrees  Celsius.  Without  cooling,  the  csmera  case 
operating  temperature  was  40  to  50  degrees  Celsius. 

With  the  optics  aligned  and  camera  case  temperature 
stabilized,  a  tare  image  (no  gas  flow  in  the  nozzle),  a 
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"flat  field"  reference  image,  and  run  laagea  were  taken  and 
atored.  For  a  run  image,  the  gaa  aupply  to  the  nozzle  waa 
turned  on  and  allowed  to  run  for  approximately  one  second. 
After  thla  start-up  period,  the  camera  control  aoftware  waa 
Initiated  caualng  the  camera  to  cycle  through  two  lmagea  at 
the  rate  of  30  Images  per  second,  then  the  gaa  aupply  waa 
turned  off.  Only  the  last  Image  of  the  two  was  retained  to 
minimize  dark  current  errors.  Thla  image  was  digitized  and 
stored  as  a  256  by  256  array  on  hard  disk  by  the  Microtex 
FORTRAN  image  management  software.  Except  for  gas  flow  in 
the  nozzle,  tare  lmagea  and  “flat  field"  reference  images 
were  taken  and  stored  similarly.  A  tare  image  consisted  of 
an  undlstorted  Airy  pattern  and  the  bias  illumination  from 
the  incandescent  light  source,  a  "flat  field"  reference 
image  consisted  of  only  the  bias  illumination  from  the 
Incandescent  light  source.  The  "flat  field"  reference 
image  was  later  subtracted  numerically  from  both  the  tare 
and  run  images.  A  set  of  tere,  reference  ,and  run  images 


waa  taken  and  etored  for  each  run  condition 


VI.  DATA 


Table  1  below.  outlines  the  jet  gas  mixture,  beam 
diameter,  and  location  for  which  Images  were  recorded.  The 
extra  Images  of  a  1  cm  diameter  beam  through  a  He  jet  at 
1  cm  downstream  were  taken  to  determine  that  the  removal  of 
one  or  more  attenuator  filters  did  not  seriously  effect 
either  the  beam  deflection  or  Strehl  ratio  measurements . 
Since  it  did  not  seem  to  greatly  effect  the  results.  this 
technique  was  used  (primarily  for  He  jets)  when  the  beam 
degradation  was  severe.  The  degradation  caused  by  SF6  jets 
however,  wss  ao  severe  that  only  a  limited  number  of  images 
using  this  gas  ware  recorded,  as  table  1  indicates.  All 
images  were  taken  with  a  1/30  second  "exposure"  time. 

The  recorded  images  may  be  visualised  using 
pseudo-three-dimensional  plots  of  the  array  response. 
Pixel  response  for  every  eighth  line,  for  example,  was 
plotted  with  the  reference  datum  for  each  line  shifted 
vertically  by  an  amount  proportional  to  the  sine  of  some 
abitrary  view  elevation  angle.  Figure  9  is  an  example  of 
such  a  plot  for  a  typical  tare  image  (no  flow)  after  the 
reference  "flat  field"  has  been  subtracted.  Ideally, 
figure  9  should  be  the  paeudo-3-D  plot  of  an  Airy 
pattern*®,  but  the  subsidiary  rings  are  not  readily 
visible.  However,  evidence  of  ring  structure  did  become 


Fig.  9.  Pixel  response  plot  of  every  eighth  row  of  s 
tare  iaage.  The  aperture  diameter  was  1  cm,  and  the 
distance  from  the  reducing  telescope  output  lens  was 
51.5  cm.  The  reference  image  of  the  bias  illumination 
has  been  subtracted  from  the  actual  data  recorded  by 


the  camera 


>  V*. 
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apparent  when  the  intensity  distribution  of  the  images  were 
Integrated  with  respect  to  radius,  as  is  shown  by  the  flat 
spots  in  the  curves  in  figures  6  and  7  and  all  subsequent 
integrations  of  intensity  for  tare  isagea.  Flat  spots 
represent  the  near  zero  contribution  to  the  integration 
froa  minima.  A  subsequent  rise  followed  by  another  flat 
spot  indicates  a  maximum.  The  integration  technique  is 
discussed  further  in  the  section  on  data  reduction  but, 
since  the  integration  is  radial,  the  maxima  and  minima  are 
inferred  to  be  ring  shaped. 


BEAN 

BEAN  DIANETER 

LOCATION 

1  cm 

0.5  cm 

He<6>,  62xHe/38xAr (3) 

He  <  3) ,  62XHe/38XAr(3) , 

x  *  1 

38XHe/62XAr (3) , 

38XHe/62xAr (3) , 

C02 C 3) ,  SFs<2> 

C02<3>,  SF& ( 3 > 

He < 3) ,  62XHe/38XAr (3) 

He ( 3) ,  62XHe/38XAr (3) , 

x  «  2 

38XHe/62XAr (3) , 

38XHe/62XAr  C  3) , 

CO2O) 

C02<3> 

Table  1.  Teat  conditions.  Numbers  in  parentheses  indicate 
the  number  of  images  recorded  for  that  combination  of  gas 
and  beam  parameters.  Beam  location,  x,  is  in  cm  downstream 


from  the  nozzle  exit 


Figures  10  and  11  were  taken  with  the  same  mas! 
aperture,  downstream  location,  at  the  same  distance  froi 
the  telescope  and  under  similar  illumination,  but  with  thi 
jet  on  and  the  beam  passing  through  the  mixing  layer 
Thus,  figures  10  and  11  may  be  compared  to  figure  9  for  ai 
appreciation  of  the  beam  degradation  caused  by  mixini 
layers  for  62*  He/38*Ar  and  CO2  jets  at  a  location  1  ci 
downstream  of  the  nozzle.  However,  each  of  the  imagei 
represented  by  figures  10  and  11  had  its  own  tare  takei 
just  previously  for  analysis  purposes. 

Of  interest  however,  is  the  fact  that  the  shape  of  thi 
plot  indicates  that  something  more  than  simple  dispersio: 
of  the  beam  is  taking  place,  the  cross  section  is  no  longe: 
circular.  As  the  distortion  caused  by  the  mixing  laye: 
becomes  more  severe,  it  also  becomes  more  complex.  Some  o 
the  images  of  degraded  beams  contained  what  appeared  to  b< 
two  peaks,  but  at  much  lower  intensities  than  thi 
undistorted  beam. 

Figure  12a  and  b  are  photographs  under  simila: 
illumination  and  with  the  same  aperture  size  and  downstreai 
location  as  the  images  plotted  in  figures  10  and  11 
Figure  12a  was  taken  with  no  flow  and  figure  12b  with  thi 
beam  passing  through  a  CO2  jet.  These  photographs  an 
provided  as  further  subjective  evidence  of  the  fidelity  o: 
the  recorded  digital  Images.  The  two  smaller  spots  in  the 
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Fig.  11.  Paeudo-3-D  Intensity  plot  for  CO2  jet,  1  cm 
diameter  beam,  1  cm  downatream  from  nozzle  exit. 


Fig*  12.  Beam  distortion  resulting  from  a  CO2  jet. 
The  center  peak  has  been  deliberately  overexposed  to 
ahow  subsidiary  ring  structure. 

12a:  Undistorted  tare. 


12b: 


1  cm  diameter  beam,  1  cm  downstream  from  nozzle 
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upper  portion  of  the  photographs  in  figure  12  are  reference 
spots  used  to  guage  the  diaplaceaent  of  the  distorted  bean 
nsxinun.  They  are  approximately  1.1  cm  apart  or  the 
original  photographs  and  are  not  themselves  displaced  from 
the  tare  to  the  jet-on  condition. 

Figures  13  through  17  are  plots  of  images  taken  with 
the  same  aize  mask  aperture,  but  at  2  cm  downstream  from 
the  nozzle.  The  pixel  response  in  figure  13  haa  been 
enhanced  by  the  removal  of  a  5x  attenuator  filter  after 
taking  its  tare,  but  before  taking  the  image  for  which 
figure  13  is  the  plot. 

Photographs,  similar  to  figure  12,  were  taken  of  the 
farfield  for  each  of  the  experimental  conditions.  While 
these  photographs  offer  little  in  the  way  oy  quantitative 
intensity  distribution  data,  they  correspond  in  overall 
structure  to  the  digital  images  that  were  recorded. 
Therefore,  the  pseudo-3-D  plota  and  the  digitized  images 
they  represent  are  believed  to  be  reasonable 
representationa  of  the  far  field  intensity  distributions  of 
the  degraded  beams. 


Fig.  13.  Paeudo-3-D  intensity  plot  for  He  jet,  1  cm 
diameter  beam,  2  cm  downstream  from  nozzle  exit.  The 
response  has  been  effectively  amplified  in  comparison 
to  other  images  at  similar  conditions  by  the  removal 
of  a  5x  attenuator  filter  prior  to  exposure. 
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Fig.  14.  Paeudo-3-D  intensity  plot  for  62XHe/38XAr 
jet,  1  cm  diameter  beam,  2  cm  downstream  from  nozzle 
exit. 


Fig.  15.  Pseudo-3-D  intensity  plot  for  38XHe/62XAr 
jet,  1  cm  diameter  beam,  2  cm  downstream  from  nozzle 
exit . 
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Fig.  16.  Paeudo-3-D  intenaity  plot  for  N2  jet,  1 
diameter  beam,  2  cm  downatraam  from  nozzle  exit. 


Fig.  17.  Paaudo-3-D  intenaity  plot  for  CO2  jet,  1  c 
diameter  beam,  2  cm  downatream  from  nozzle  exit. 
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dimensions  and  because  of  the  odd-even  response  of 
alternate  lines,  but  becauaa  it  waa  found  that  aoae  aspect 
of  the  image  capture  and  digitization  system  caused  a  drift 
in  the  zero  datum  level  of  successive  frames  that  could  not 
be  quantified  on  a  pixel  by  pixel  baaia.  In  order  to 
correct  for  this,  the  following  approach  waa  used: 

1.  The  location  of  the  intensity  centroid  of 
the  pixels  having  recorded  intensity  values  not  less  than 
90*  of  the  maximum  value  recorded  for  the  image  was 
determined  numerically. 

2.  The  curves  like  those  plotted  in  figures  6 
and  7  were  developed  by  a  summation  process  of  the  form. 


P(mAr)  =221. 

s  t  Sl 


-  s  )*  +  ( 1  -  t  )2  <  m  Ar 


<9> 


where: 

Z«t  *  the  recorded  intensity  at  the  pixel  located  in 
the  a'th  row  and  the  t ' th  column. 

Ikl  =  the  intensity  of  the  pixel  located  nearest  the 
of  the  centroid  of  the  peak  of  the  intensity 


distribution 
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k  =  the  row  number  of  the  location  of  the  centroid  of 
the  peak  of  the  intenaity  diatribution . 

1  *  the  column  number  of  the  location  of  the  centroid 
of  the  peak  of  the  intenaity  diatribution. 

*  *  1,  2,  3,  ...»  1000. 

At  *  aindjaaxr  3 max) /1000. 

iiiax  *  leaat  number  of  rowa  between  k  and  either  edge 
of  the  array. 

Jmax  =  leaat  number  of  columna  between  1  and  either 
the  top  or  bottom  of  the  array. 


f 


Thia  summation  approximated  the  Integral  relationship. 


l 

H 


P(r)  = 

0 


i 

«| 


Since 


Ke)  2  7T  e 

i 


where : 


P  =  is  the  value  of  P  for  the  degraded  bean, 
Pt  *  the  value  of  P  for  the  tare. 


3.  The  Unit  in  eqn.  <11>  was  then  estimated 
graphically  by  plotting  the  ratio 


P(nAx) 

Pt(nAx) 


veraua  the  nondimensional  radius,  x,  given  by: 


x  =  kaw 


where: 


k  *  the  wave  number, 
a  =  the  diameter  of  the  beam  waist. 


xmmaam 


JUHIUIIIW 
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w  =  r/D,  the  ratio  between  the  radius  and  the 
distance  fro*  the  bee*  waist. 

It  was  important  to  estimate  the  Strehl  ratio 
graphically  because  of  the  extreme  variations  in  the  ratio 
in  the  limit,  eqn.  (11),  as  r  approached  Ar.  These 
variations  are  due  to  the  non-continuous  nature  of  the 
digitized  data  as  a  function  of  the  non-continuous 
parameter  nr.  As  an  example  of  how  this  was  done,  two 
separate  images  of  the  same  undegraded  beam  were  recorded 
and  the  Integrated  intensities  and  ratio  of  Integrated 
intensities  are  plotted  in  figure  18. 

Figure  18  shows  three  curves.  The  two  curves  which 
are  initially  lowest  are  plots  of  the  normalized,  numerical 
Integration  of  the  beam  intensities  over  the  radius  as 
functions  of  radius.  The  upper,  nearly  horizontal,  curve 
is  the  ratio  of  the  the  integrated  intensities  as  a 
function  of  radius. 

In  this  figure,  the  Strehl  ratio  of  two  images  taken 
from  the  same  undegraded  beam,  but  at  different  times,  is 
estimated  from  the  plot  of  the  ratio  of  the  integrated 
intensities,  or  powers,  of  each  beam  as  a  function  of 
radius.  This  curve  is  the  numerical  approximation  of  the 
ratio  in  expression  (12).  The  theoretical  value  of  the 
Strehl  ratio  in  this  case  la,  of  course,  1.0,  and  the  graph 


t  aw  n/  TPTT  ~  CT  i  T't*  -  J  :  V.l  FS'f  ‘g  ,  W  l»J  ■■  f  a^JI  "JJ  V»-C»  in  7TW fj  LTIW.W.T  ^ 
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can  be  seen  to  approach  that  value  in  apite  of  the  "noise" 
near  the  Halt  ee  r  approachee  0. 

Another  example  of  thia  graphical  technique  for 
estimating  Strehl  ratio,  uaing  actual  degraded  beam  data, 
ia  ahown  in  figure  19.  For  reference,  thia  ia  the  plot 
uaed  to  estimate  the  Strehl  ratio  for  the  iaage  whoae 
paeudo-3-D  intenaity  plot  ia  ahown  in  figure  10.  A  similar 
plot  waa  generated  for  each  digitized  image  that  was 
recorded.  The  Strehl  ratios  for  the  recorded  images  were 
obtained  by  estimating  the  value  that  the  plot  of  the  power 
ratio  tended  toward  aa  the  value  of  the  radiua  became 
small . 

Tables  2  and  3  present  a  summary  of  the  results  of  the 
analysis  of  all  of  the  recorded  data.  For  each  combination 
of  experimental  conditions,  the  average  of  the  estimated 
Strehl  ratios  ia  tabulated. 


Fig.  19.  Normalized  power  and  power  ratio  plots  for 
an  undegraded  tare  and  a  1  ca  beam  1  ca  downatreaa  in 
a  62*He/36*  He  jet.  This  is  the  radial  integration  of 
figure  10. 


.Power  ■}  Power  ratio  vs  nondi^ensionel  rsoius 


APERTURE  DIAMETER 


0.5  cm 

1  cm 

DENSITY 

AVG  STREHL 

RATIO 

AVG  STREHL 

RATIO 

RATIO 

7.22 

0.22 

0.07 

1.64 

0.40 

0.15 

1.10 

0.66 

0.34 

N2 

1.03 

1.0 

0.97 

C02 

0.65 

0.44 

0.20 

VIII.  DISCUSSION 


The  estimated  Strehl  ratio  values  from  the  data  are 
plotted  versus  An2c2x2  figUre  20  along  with  curves  for 
eqn.  <6>  using  values  for  a  =  0.1S  and  a  =  0.20.  The 
value  for  c  used  for  eqn.  <6>  in  these  plots,  as  a  function 
of  density  ratio  was  determined  from  the  three  data  points 
on  the  spreading  rate  given  by  Brown  and  Roshko*0  and  the 
following  curve  fit  function: 

Pi  Pi 

c  *0.38 +  0.14  log - +log2 -  (14> 

Pi  Pi 

Note  that  eqns.  <5>  and  <6>  predict  that  aperture  size 
should  have  no  effect  on  the  Strehl  ratio  since  aperture 
size  does  not  appear  in  the  equation.  However,  examination 
of  figure  20  reveals  that  for  each  value  of  An2c2x2,  the 
measured  data  tends  to  fall  into  two  groups:  a  higher 
value  of  Strehl  ratio,  and  a  lower  value.  The  higher 
grouping  always  corresponds  to  the  smaller  aperture  size, 
and  the  lower  grouping  to  the  larger  aperture. 
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It  should  also  be  noted  that  the  approach  used  to 
estinate  the  Strehl  ratio  assumes  that  the  peak  Intensity 
centroid  la  the  center  of  the  far  field  focus,  aa  does  the 
measurement  of  beam  deflection  mentioned  below.  This  may 

not  be  the  case,  but  should  be  a  cloae  approximation 
to  the  tilt  correction  required  by  the  definition  of  the 
Strehl  ratio.  Furthermore,  in  the  cases  where  the  beam  was 
no  longer  even  circular  its  not  clear  how  one  would  define 
the  center  of  the  far  field  focus. 

Measured  beam  deflection  data  for  all  of  the  test 
conditions  are  plotted  in  figure  21.  Negative  values  of  An 
(DELTA  n  on  the  plot)  indicate  jet  gaaea  with  indices  of 
refraction  less  than  that  of  air,  while  positive  values 
indicate  that  the  jet  had  a  larger  index  of  refraction  than 
air.  The  deflection  angle  is  the  angle  the  beam  waa  tilted 
away  from  the  nozzle,  aa  determined  by  the  movement  of  the 
location  of  peak  intensity.  Hence,  negative  deflection 
angles  indicate  that  the  beam  was  deflected  toward  the 
nozzle . 

It  is  interesting  to  note  that  eqn.  (4)  predicts  that 
a  homogeneous  wedge  perpendicular  to  the  beam  will  not 
deflect  the  beam.  But,  as  the  wedge  is  tilted,  eqn.  (4) 
predicts  increasing  beam  deflection.  Further,  if  the  tilt 


of  the  wedge  is  away  from  the  jet  at  an  angle  of 
approximately  20*  of  the  total  wedge  angle  the  beam 
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deflections  predicted  by  eqn .  (4)  are  very  close  to  the 
Measured  bean  deflections.  Such  a  tilt  is  consistent  with 
the  shapes  of  mixing  layers  seen  by  Poling^  using  the  same 
nozzle  as  wss  used  here. 


IX.  GQfcJCkUSIQNS 


The  Strehl  ratios  measured  in  this  work,  and  plotted 
in  figure  20,  show  roughly  the  same  variation  with  An^c^x^ 
as  is  predicted  by  Sutton  in  eqn.  (6),  except  for  the  He 
jet  data  taken  at  2  cm  downstream.  It  is  likely  that  a 
value  of  a  could  be  found  for  each  of  the  other  test 
conditions  which  would  make  the  match  for  the  remaining 
data  even  closer,  as  it  is  possible  that  the  turbulence 
intensity  varied  from  jet  to  jet.  However,  a  was  not 
measured.  And,  it  is  generally  held3  that  the  presence  of 
coherent  structure,  as  was  the  case  here,  reduces  the 
validity  of  eqn.  (6)  as  a  predictor  of  beam  degradation. 
Further,  the  data  demonstrate  a  variation  with  beam 
diameter  not  predicted  by  eqn.  (5>  or  eqn.  (6) . 

The  He  jet  data  considered  to  be  the  exception  are  the 
data  at  the  far  right  in  figure  20.  The  intensity 
distribution  for  the  beam  after  passing  through  the  He  jet 
2  cm  downstream  is  nearly  indistinguishable  from  the 
background,  even  after  the  removal  of  a  5x  attenuator 
filter.  This  is  shown  by  figure  13  which  is  the  pseudo-3-D 
intensity  plot  for  one  of  these  images.  These  data  are 
considered  to  be  highly  suspect  because  the  recorded  data 
is  so  close  to  the  noise  level  of  the  array  output. 
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The  measured  beam  deflections  also  show  a  trend 
similar  to  a  prediction  model,  that  of  a  simple  wedge  with 
two  optical  interfaces.  However  the  wedge  in  this  simple 
model  muat  be  tilted  away  from  the  jet  and  toward  the 
Incoming  beam  in  order  to  match  the  data.  This  seema  to  be 
a  reasonable  requirement  baaed  on  the  observed  shape  of 
actual  mixing  layers. 

Here,  an  arbitrary  tilt  angle  for  the  wedge,  equal  to 
20*  of  the  total  wedge  angle,  was  chosen  because  it  seems 
natural  that  the  amount  of  tilt  would  depend  on  the  same 
parameters  that  determine  the  wedge  angle,  and  the  factor 
of  20*  fit  the  data.  more  complex  model  using  a  finite 
number  of  small  wedges  may  do  a  better  job  of  modeling  the 
beam  deflection.  This  would  be  a  problem  of  quadrature 
though,  because  the  beam  angle  of  Incidence  for  each  of  the 
wedges  would  depend  on  the  angle  of  transmission  from  the 
preceding  wedge. 

In  all,  the  results  of  this  work  validate  the  concept 
of  photo  sensor  arrays  as  sensors  to  directly  sample  the 
far  field  distortion  caused  by  turbulent  interfaces.  And, 
it  provides  direct  experimental  data  verifying  in  general 
the  dependencies  identified  by  Sutton17  and  Vu  et.  al6. 
However,  it  suggests  the  possibility  of  other  dependencies 
associated  with  the  size  of  the  beam  relative  to  soma 


length  dlmenaion  of  the  mixing  layer.  The  nature  of  this 
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dependency  wae  not  discovered. 

Further  work  on  the  development  of  this  technique  and 
it's  use  to  expand  the  range  of  experimental  conditions, 
primarily  various  density  ratios  at  higher  Mach  numbers,  is 
warranted . 
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